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ABSTRACT 



Aims. The interstellar medium (ISM) has a multiphase structure characterized by gas, dust and molecules. The gas can be found in 
different charge states: neutral, low-ionized (warm) and high-ionized (hot). It is possible to probe the multiphase ISM through the 
observation of its absorption lines and edges in the X-ray spectra of background sources. 

Methods. We present a high-quality RGS spectrum of the low-mass X-ray binary GS 1826-238 with an unprecedent detailed treat- 
ment of the absorption features due to the dust and both the neutral and ionized gas of the ISM. We constrain the column density 
ratios within the different phases of the ISM and measure the abundances of elements such as O, Ne, Fe and Mg. 
Results. We found significant deviations from the proto-Solar abundances: oxygen is over-abundant by a factor 1.23 ± 0.05, neon 
1.75 + 0.11, iron 1.37 ±0.17 and magnesium 2.45 ±0.35. The abundances are consistent with the measured metallicity gradient in our 
Galaxy: the ISM appears to be metal-rich in the inner regions. The spectrum also shows the presence of warm/hot ionized gas. The 
gas column has a total ionization degree less than 10%. We also show that dust plays an important role as expected from the position 
of GS 1826-238: most iron appears to be bound in dust grains, while 10-40% of oxygen consists of a mixture of dust and molecules. 
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The interstellar medium of our Galaxy (ISM) is a mixture of 
dust and gas in the form of atoms, molecules, ions and elec- 
trons. It manifests itself primarily through obscuration, redden- 
ing and polarization of starlight and the formation of absorption 
lines in stellar spectra, and secondly through various emission 
mechanisms (broadband continuum and line emission). The gas 
is found in both neutral and ionized phases (for a review, see 
Ferriere 2001). The neutral phase is a blend of cold molecular 
gas (T ~ 20-50 K), found in the so called dark clouds, and cold 
atomic gas (T ~ 100 K) inherent in the diffuse clouds, while 
the warm atomic gas has temperatures up to 10 4 K. The atomic 
gas is well traced by H i and mainly concentrated in the Galactic 
plane with clouds up to few hundreds pc above it. The warm ion- 
ized gas is a low-ionized gas, with a temperature of ~ 10 4 K. It 
is mainly traced by Ha line emission and pulsar dispersion mea- 
sures; it can reach a vertical height of 1 kpc. The hot ionized gas 
is characterized by temperatures of about 10 6 K. It is heated by 
supernovae and stellar winds from massive stars; it gives rise to 
high-ionization absorption lines and the soft X-ray background 
emission. The study of the ISM is very interesting because of its 
connection with the evolution of the entire Galaxy: stellar evolu- 
tion enriches the interstellar medium with heavy elements, while 
the ISM acts as source of matter for the star forming regions. 

High resolution X-ray spectroscopy has become a powerful 
diagnostic tool for constraining the chemical and physical prop- 
erties of the ISM. Through the study of the X-ray absorption 
lines in the spectra of background sources it is possible to probe 
the various phases of the ISM of the Galaxy. First of all, the K- 



shell transitions of low-Z elements, such as oxygen and neon, 
and the L-shell transitions of iron fall inside the soft X-ray en- 
ergy band. Secondly, the presence of different charge states for 
each element allows us to constrain the multiphase ISM, e.g. its 
ionization state and temperature distribution. 

Schattenburg & Canizaresl ([1986) first measured ISM ab- 
sorption edges in the X-ray band using the Einstein Observatory 
and they found features consistent with the Oi Is - 2p line 
and traces of On. After the launch of the XMM-Newton and 
Chandra satellites a new era for the study of the ISM opened up. 
The grating spectrometers on board of these satellites, RGS and 
LETGS/HETGS respectively, provide a spectral resolution high 
enough to resolve the main absorption edges and lines. Recently, 
Yao et al. (2009) found high-ionization absorption lines of ions 
such as O vi to O vra and Ne vm to Ne x in the HETGS spec- 
trum of the low-mass X-ray binary Cyg X-2, and argued that 
the bulk of the O vi should originate from the conductive in- 
terface between the cool and the hot gas. Other work has re- 
vealed complex structure around the oxygen K-shell absorption 
edge (Paerels et al. 200li: lde Vries et aLll2003t I Juett et al.ll2004l) . 
ICostantini et al.l d2005l) argued that the feature of the scattering 
halo of Cyg X-2 near the O i K-edge can be attributed to the 
presence of dust towards the source, with a major contribution 
from silicates such as olivine and p yroxene. In their paper on Sco 
X- 1 , observed with XMM-Newton. \de Vries & Costantinil (2009) 
found clear indications of extended X-ray absorption fine struc- 
tures (EXAFS) near the absorption edge of oxygen. 

In this work we report the detection of absorption lines 
and edges in the high-quality spectrum of the low-mass X-ray 
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Table 1. Observations used in this paper. We report the total ex- 
posure length together with the net exposure time remaining af- 
ter screening of the background and removal of bursts. 



ID 



Date 



Length 
(ks) 



RGS 

(ks) 



PN 
(ks) 



0150390101 
0150390301 



2003 April 6 
2003 April 8 



108 

92 



67.8 
77.8 



63.8 
67.8 



binary (LMXB) GS 1826-238 obtain ed by the XMM-Newto n 
Reflection Grating Spectrometer (RGS, Iden Herder et al.ll200lh . 
In order to const rain the continuum parameters we also 
used the EPIC-p n (Striideret al. 2001) dataset of this source. 
iThompson ela l. (2008) using XMM-Newton and RXTE obser- 
vations on April 2003 derived a high unabsorbed bolometric flux 
F ~ 3.5 x 10~ 12 W rrT 2 . The source is well suited for the analysis 
of the ISM also because of its column density Nh ~ 4 x 1 25 irT 2 
(see Table [3]), high enough to produce prominent O and Fe 
edges. We assume the distance of the source to be 6.1 + 0.2 kpc 
dHeger et all2007l) . 

We analyze the absorption in the spectrum as follows. We 
first remove the bursts, because they add a strongly variable com- 
ponent to the spectrum. Then we determine the source contin- 
uum by fitting simultaneously EPIC and RGS data. In second in- 
stance we use only the high-resolution RGS spectra to constrain 
the absorption contributions. We search for statistically signifi- 
cant features by adding several absorbers in sequence: cold gas, 
warm gas, hot gas, dust and molecules. All of these appear to be 
important. 



2. Observations and data reduction 

GS 1826-238 (Galactic coordinates I = 9°.27, b = -6° .09) has 
been observed twice with XMM-Newton for a total length of 200 
ks (see Table[TJfor details). The data are reduced with the XMM- 
Newton Science Analysis System (SAS) version 9.0.1. 

GS 1826-238 is a bursting LMXB with a regular time sepa- 
ration between the bursts. Because the primary aim of the XMM- 
Newton observations was the study of the bursts, the EPIC-pn 
detector was operated in timing mode, which means that imag- 
ing is made only in one dimension, along the RAWX axis. Along 
the row direction (RAWY axis), data from a predefined area on 
one CCD chip are collapsed into a one-dimensional row for a 
fast read-out. Then source photons are extracted between RAWX 
values 30 - 45 and background photons are extracted between 
rows 2 - 16, as recommended by the standard procedure. 

We produce pn lightcurves mainly to remove the bursts in- 
tervals and to extract the spectra of the persistent part of the 
lightcurve. In the first observation 9 bursts have been detected, 
in the second observation 7 bursts. We plot the burst profiles of 
these 16 bursts in Fig. [2] We estimate a mean duration of about 
300 s for the bursts, and we remov e for each bu r st 50 s before the 
peak to 250 s after it. Recentlv lin't Zand et al.1 (|2P09) suggested 
a mean duration of about 1 ks for the bursts, but they also argued 
that after the first 100 s the inferred emission decreases sharply 
by at least an order of magnitude, contributing only about 3% to 
the fluence in the burst. After 250 s the flux of the burst has de- 
creased by almost 2 orders of magnitude and its profile merges 
with the persistent lightcurve. Thus, by removing 300 s for each 
burst, we retain less than ~ 1% burst emission, which is negligi- 
ble compared to the persistent emission. 

We process RGS data with the SAS task rgsproc. We pro- 
duce the lightcurves for the background in CCD9 following the 
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Fig. 1. Lightcurves of the first (left) and the second (right) obser- 
vation with RGS 1 . Intervals with high background have already 
been taken out. The bursts are still shown for displaying purpose, 
but their contribution in the spectra are removed as described in 
the text. The plot shows the quasi-periodicity of the bursts and il- 
lustrates why GS 182 6-238 is called the 'clock-burster' LMXB 
(lUbertini et alJl9 99). 
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2. Mean profile of the bursts in the RGS lightcurve of the 
observation. The zero point of the time scale is centered on 
burst profile peak. The red line represents the mean count 
of the 2 ks around the peaks. 



XMM-SAS guideQin order to remove soft proton flares and spu- 
rious events. We create good time intervals (GTI) by removing 
intervals with count rates higher than 0.5 c s . We reprocess 
the data again with rgsproc by filtering them with the GTI for 
background screening and bursts removal. We extract response 
matrices and spectra for the two observations. The final net ex- 
posure times are reported in Table [TJ 

Our analysis focuses on the 7 - 3 1 A (0.4 - 1 .77 keV) first 
order spectra of the RGS detector. In order to fit the spectral 
continuum properly, we also use the 0.5 - 10 keV EPIC spectra 
of both observat ions. We perform s pectral analysis with SPEX0 
version 2.01.05 (Kaastr a et al.lll996l) . We scale element al abun- 
dances to the proto-Solar abundances of Lodders (2003): N/H = 
7.943xl0~ 5 , O/H = 5.754xl0~ 4 , Ne/H = 8.912xl0 -5 , Mg/H = 
4. 169xl0~ 5 , Fe/H = 3.467xl0~ 5 . We use the C-statistic through- 
out the paper and adopt \cr errors. 



1 http://heasarc.nasa.gov/docs/xmm/abc/ 

2 www.sron.nl/spex 
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3. Spectral Modeling 

The first step of the spectral analysis consists of the determi- 
nation of the continuum emission and the dominant absorp- 
tion component. The best way to do this is to fit the spectra 
of RGS and EPIC-pn simultaneously. The XMM-Newton cross- 
calibration is very complex, not only because of the different en- 
ergy bands, but mainly because of their different features: RGS is 
sensitive in the soft X-ray energies with high spectral resolution, 
showing narrow absorption features, while pn has a low spectral 
resolution, therefore blurring the absorption features seen with 
RGS. EPIC-pn has a higher count rate compared to RGS and a 
broader energy band. The original RGS spectra are binned by a 
factor of 10 in this simultaneous fit. This is necessary to tem- 
porarily remove the narrow features due to the absorption lines. 
The pn spectra are resampled in bins of about 1/3 of the spectral 
resolution (FWHM ~ 50-150 eV between 0.5-10 keV), which 
is the optimal binning for most spectra. 

A better local fit for absorption edges and lines is obtained 
from a separate RGS fit. In the RGS local fit we rebin the spectra 
only by a factor of two, i.e. about 1/3 FWHM (the first order 
RGS spectra provide a resolution of 0.06-0.07 A). This gives at 
least 10 counts/bin and a bin size of about 0.02 A. 

3.1. Simultaneous EPIC-RGS fits 



Table 2. EPIC-RGS spectral fits to the persistent emission. 



At firs t we follow the spectral modelling of Thompson et al.l 
(2008). The continuum spectrum is modeled by emission from 
a black body and two comptonization models. The black body 
component arises from the thermal emission of the accretion 
disk around the neutron star. The first comptonization compo- 
nent (hereafter CI) describes the energy gain of the disk soft 
photons by scattering in the accretion disk corona. The second 
comptonization component C2 corresponds to scattered seed 
photons originating from regions clo s er to the NS surface, i.e. 
the boundary layer. Thomp son et alj (I2008I) applied a neutral 
absorber to the continuum mentioned above and fitted XMM- 
Newton, Chandra and RXTE data. For this purpose we use the 
absm model in SPEX: the model calculates the continuum trans- 
mission of neutral gas wi th cosmic abundances as published by 
iMorrison & McCa mmonl d!983l) . In our case the same model 
does not give a satisfactory fit, especially around the neon and 
oxygen edges, and cannot fit the Oi line. This could be ex- 
pected because the Morrison & McCammon ( 1983) model does 
not take into account the absorption lines and the possible vari- 
ations in the abundances. Therefore we replace the absm com- 
ponent with a hot component, which describes the transmission 
through a layer of collisionally ionized plasma. At low temper- 
atures it calculates the absorption of (almost) neutral gas (for 
further informations see the SPEX manual). We leave the tem- 
perature and the O, Ne, Mg and Fe abundances of this absorber 
free in the fit. In the fits we ignore two small regions (17.2-17.7 
A and 22.7 - 23.2 A), close to the iron and oxygen edges respec- 
tively. The presence of dust and molecules affects the fine struc- 
ture of the edge, thus these regions will be analyzed with more 
complex models in Sect. 13.2.31 However the ISM abundances 
are determined by the depth of the absorption edges, thus ig- 
noring these small regions we still can constrain the abundances 
of such elements dKaastra et al.ll20 09). Indeed, in Sect. I3.4l and 
Table [6] we will validate this assumption. We obtained a good fit 
with C-stat/dof0 = 2451/1705 and 2579/1710 in the two ob- 
servations (see Fig. [3}. The parameters for both the observations 



Par / component 


OBS 1 


OBS 2 


Average value 


absorber 








N H (10 25 m- 2 ) 


3.65 ± 0.05 


3.68 ± 0.04 


3.67 ± 0.03 


kT (l(T 4 keV) 


7.07+0.10 


7.22 ± 0.09 


7.15 ±0.07 


O 


1.471 ±0.008 


1.403 ±0.008 


1.437 ±0.006 


Ne 


2.72 ± 0.04 


2.64 ± 0.04 


2.68 ± 0.03 


Mg 


0.81 ±0.09 


0.80 ± 0.09 


0.80 ± 0.06 


Fe 


1.90 ±0.03 


2.01 ±0.02 


1.98 ± 0.02 


black body 








nux(l(T 13 WrrT 2 ) 


0.65 ±0.10 


0.55 ± 0.03 


0.56 ± 0.03 


kT bb (keV) 


0.170 ±0.002 


0.167 ±0.002 


0.168 ±0.001 


CI comptonization 








flux(10- 13 Wrrr 2 ) 


7.8 ± 0.3 


7.7 ± 0.3 


7.75 ± 0.21 


kT, (keV) 


0.25 ± 0.02 


0.24 ± 0.01 


0.24 ± 0.01 


kT s (keV) 


2.07 ± 0.05 


2.38 ± 0.06 


2.20 ± 0.04 


T 


11.0 + 0.6 


9.3 ± 0.4 


9.8 ±0.3 


C2 comptonization 








flux(l(r 13 Wrrr 2 ) 


0.9 ±0.1 


0.6 ±0.1 


0.75 ± 0.07 


kT, (keV) 


0.50 ± 0.01 


0.54 ± 0.01 


0.52 ±0.01 


kT c (keV) 


9.4 ±0.8 


4.6 ± 0.5 


5.9 ±0.4 


T 


<0.5 


<0.7 


- 


C sm 1 dof 


2451/1705 


2579/1710 


- 



Here and hereafter dof means degrees of freedom. 



Abundances are relative to the proto-Solar values of iLoddera (2003). 
Fluxes are derived in the 0.3-10 keV band. We also report the weighted 
averages between the two observations. See also Fig. [3] 



are listed in Table [2] We designate this simple model where the 
ISM is modeled with one (neutral) gas component as Model A. 
The abundances are mostly in agreement between the two ob- 
servations, but they are not trustworthy. In Sect. 13.21 and !3. 41 we 
show that the RGS fit provides a column density higher by 10%, 
which significantly changes the abundances estimates. There are 
also small differences in the continuum parameters, such as the 
electron temperatures. Indeed we find different temperatures for 
both comptonization components (see Table|2]). These small de- 
viations affect the broadband spectral slope and forbids to fit the 
two EPIC-pn observations simultaneously, while this is possible 
with the RGS spectra. 

We also test alternative continuum models in order to show 
that the adopted model is the best one. iThompson et all ([2008) 
showed that the spectral modeling of GS 1826-238 can be done 
with other continuum models: 1) black body emission plus a sin- 
gle comptonization component; 2) black body emission plus a 
cut-off powerlaw; 3) double comptonization plus a disk black 
body; 4) two comptonization components. We test them on both 
spectra, but report here only the results for the first spectrum. 
The results for the other observation are similar. The models (1) 
and (2) give similar results, but with C-stat/dof = 3000/1709 the 
fit is worse than for our adopted Model A. Model (3) gives even 
worse fits. The final alternative model (4) gives an intermediate 
result C-stat/dof = 2625/1707. 

3.2. The high-resolution RGS spectra 

In Fig. [4] we plot the RGS spectrum of the persistent emission 
in the first and second observation. Several interesting features 
can be recognized. At 23.1 A we see the absorption edge of the 
interstellar neutral oxygen, while the Oi line is clearly visible 
at 23.5 A. There is also a broad absorption feature close to the 
oxygen edge, which is clearly seen in the fits residuals, see Sect. 
I3.2.3l for a dedicated discussion. The K-edge of neon and L-edge 
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Fig. 3. Simultaneous RGS-EPIC best fits of the first (top) and 
the second (bottom) observations. The model used is Model A 
(see Sect. I3.ll and Table [2}. The upper panels show from top 
to bottom the EPIC and RGS spectrum, respectively. The lower 
panels show the fit residuals (dark points with small error bars: 
EPIC). The dips present in the count spectrum of the RGS cor- 
respond to bad columns in the RGS with lower sensitivity. 



of iron are easily recognized at 14.3 A and 17.5 A, respectively. 
We fit the RGS spectra of the two observations simultaneously 
with Model A, freezing the shape of the continuum emission and 
leaving as free parameters the normalizations of the emission 
components and the parameters of the absorber. In the fits we 
still ignore the two small regions (17.2 - 17.7 A and 22.7 - 23.2 
A), close to the iron and oxygen edges respectively (see Sect. 
I3.2.3l for the dedicated analysis). The results of the RGS spec- 
tral fits are shown in Tableland are designated as Model A. We 
report also the results of the fits obtained for each observation: 
the agreement between the parameters validates the simultane- 
ous spectral fit. As expected, the residuals (Fig. |4]l show large 
deviations in the spectral regions that we have temporarily re- 
moved: ~ 4<x and ~ 3<x near 17.4 A and 23 A, respectively. 
These features cannot be modeled with a pure-gas model and re- 
quire the introduction of dust and molecular components in our 
model (see Sect. 13.2.31 . 



3.2.1. The neutral gas 

The fits obtained with a simple model (a single gas component 
for the ISM) show that the ISM can be initially modeled with 
cold gas (see Table [2] and [3] Model A). It has a mean temper- 
ature of kT ~ 6 x 1CT 4 keV, i.e. about 7 000 K, and provides 
the bulk of the warm atomic gas anticipated in the introduction. 
It is almost neutral, except for iron and magnesium: Mg n con- 
tributes 30 % to the total magnesium column density, while Fe n 
accounts for 20 % of the iron. A precise measure of the ratios 
Mg i/Mg n and Fe i/Fe n for our spectra is not possible. Near the 
magnesium edge the spectrum is noisy and the Mg i and Mg n 
edges are close, at 9.48 A and 9.30 A respectively, while near 
the iron edge the lines are unresolved and there is also a contri- 
bution from dust that affects the edge structure (see Sect. l3.2.3l >. 
However, the total magnesium and iron column densities are es- 
timated taking into account the jump across the respective edges 
and they will not be affected by these problems. As expected, the 
RGS spectral fits provide a different A^h value than the simul- 
taneous EPIC-RGS fit, because of the imperfect EPIC-RGS 
cross-calibration and the low resolution of EPIC, that smooths 
the absorption features (see Fig . [3j ■ 



3.2.2. The ionized gas 

The fit residuals near 21.6 A and 23.35 A (see Fig. @}, where 
we should expect the ls-2p transitions of Ovn and On respec- 
tively, suggest the presence of additional weak absorption fea- 
tures. Other weak features are found near 13.4 A and 14.6 A, 
close to the theoretical Ne ix and Ne n wavelengths. We deal sep- 
arately with the different ionization states. 

At first we make a fit to the RGS spectra adding columns 
of On and Nen to our model through a slab component. The 
slab model calculates the transmission of a layer of plasma with 
arbitrary composition. Free parameters are the intrinsic veloc- 
ity dispersion and t he column densities of the individual ions 
dKaastra et al.lll996h . The fits improve significantly: by fitting 
the two RGS observations simultaneously we get AC-stat ~ 130. 
The velocity dispersion is not well constrained (cry — 50 + 15 
km s _1 ). The average ion columns are 1.2 + 0.4 xlO 21 irT 2 (On) 
and 2.4 ± 0.4 xlO 21 irT 2 (Ne n), while the cold gas gives 3.05 + 



0.15 xl0 22 irr z (Oi) and 6.7 + 0.3 xlO 21 mT l (Nei). 

In second instance we add another slab component to take 
into account the contribution by the hot gas . The columns are 1 . 1 



19 ™ -2 



+ 0.5 xlO^m^ (both O vn and O vm) and 3.5 + 2.5 xl0 ly m 
(Ne ix). The addition of the hot gas provides AC-stat = 30, which 
is significantly less than the improvement we have obtained by 
adding the low-ionized gas. Moreover we can only put an upper 
limit to the velocity dispersion of the hot ionized gas (250 km 

However, in order to take care of every absorption feature 
created by all ions in the warm-hot phases and to deal with 
physical models, we substitute the two slab components with 
two hot components. We couple the elemental abundances of 
the warm-hot components to those of the cold gas, assuming all 
ISM phases have the same abundances. This is a reasonable as- 
sumption, especially for the warm (low-ionization) ionized gas, 
as its temperature is not too different from the temperature of 
warm neutral gas. The additional free parameters are the hydro- 
gen column density, the temperature and the velocity dispersion. 
In summary, the additional warm and hot phases give an average 
improvement of AC-stat ~ 80 for only 6 free parameters added. 
We label such a 3-gas model as Model B and display all results 
in Table [3] We plot the individual absorption edges of O, Fe, Ne 
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Fig. 4. Continuum best fit to the RGS spectrum of the first (black) and the second (orange) observation. Here we use the simple 
model used for EPIC-RGS fits (see Sect. 13. 11 1. In the fits we excluded two small regions near the Oi K-edge and the Fei L-edge 
(see Sect. l3.2l ). which are indicated by two red horizontal strips in the top panel. See Sect. l3.2.3l for the detailed analysis. The results 
of the fits are shown in Table [3] they refer to model A. 



and Mg in Fig. |5J |6] [7] and [U respectively. We discuss these re- 
sults in Sect. 13.31 The predicted deviations near the oxygen and 
iron edges, clearly seen in Fig. [5] and [6l still confirm that pure 
interstellar gas cannot reproduce all the absorption features and 
that we need to take into account different states of matter, such 
as solids. 

Finally, we also consider the model where the cold gas com- 
ponent is forced to be neutral by freezing its temperature to 
5 x 1CT 4 keV, i.e. 5 800 K. In this case, we get an almost equal 
fit (C-stat/dof = 4466/3227), with the exception that the warm 
component has a significantly lower temperature (~ 1.4 x 10~ 3 
keV, i.e. 10-20 000 K), to account for the O n that in our nom- 
inal fit is partially produced by the cold component. This tem- 
perature value is more representative than our previous va lue of 
5.4 x 10~ 3 keV for the warm ionized gas in the ISM (iFerrierd 
2001). However, both the fits are acceptable, thus we report only 
results obtained with the cold-gas temperature as free parameter 
in Tabled 

3.2.3. Fine structures: Dust and molecules 

Further important improvements to our fit are obtained by adding 
both dust and molecules to our multiphase gas model. We use 
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Fig. 5. 1 K-edge: data and Model B. Black and grey points refer 
to the first and second observation, respectively. 

here two recently added models of SPEX: dabs and amol, which 
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Table 3. RGS spectral fits to the persistent emission. 







Mod A (fl) 


ModB 


ModC 


Component 


Parameter 


OBS 1 


OBS 2 


OBS 1+2 


OBS 1+2 


OBS 1+2 


Cold 


A^OO^m- 2 ) 

kT (10~ 4 keV) 

o> (km s _1 ) 

O 

Ne 

Mg 

Fe 


4.18 ±0.01 
5.90 ±0.14 
27 ±15 
1.30 ±0.02 
2.08 ± 0.04 
2.27 ±0.12 
1.39 ±0.02 


4.22 ± 0.07 

6.13 ±0.12 
18 ±6 
1.29 ±0.02 
1.86 ±0.07 

2.14 ±0.16 
1.42 ±0.06 


4.21 ± 0.09 
6.02 ± 0.09 
13 ±7 
1.29 ± 0.01 
1.95 ±0.07 
2.21 ±0.16 
1.39 ±0.05 


3.88 ± 0.07 
6.04 ±0.10 
< 12.6 
1.29 ±0.02 
2.19 ±0.10 
1.93 ±0.15 
1.65 ±0.08 


3.94 ± 0.05 

8.6 ± 0.4 

<24 

1.17 ±0.03 

1.75 ±0.11 

1.30 ±0.25 

<0.1 


Warm 


^H(10 i5 m- 2 ) 
kT (10~ 3 keV) 
o> (km s _1 ) 


- 


- 


- 


0.46 ± 0.06 
5.4 ±0.3 
50 ±25 


0.15 ±0.05 
4.5 ± 0.5 
< 150 


Hot 


A^OO^m" 2 ) 
kT (keV) 
o"v (km s ) 


- 


- 


- 


0.042 ± 0.008 
0.20 ± 0.02 
< 160 


0.047 ±0.011 
0.20 ± 0.03 
<200 


Dabs (i) 


N Fe 


" 


" 


~ 


" 


2.6 ±0.1 
2.3 ±0.1 


Amol (Ac) 


N (i=14, Silicates) 
N (i=7, H 2 Ice) 
No (i=2, CO) 
No (i=23, Aluminates) 


- 


- 


- 


- 


2.5 ±0.5 
<0.7 
<0.4 
<0.4 


Statistics 


C-Stat / dof 
C-Stat / dof (*) 


2239/1595 
2930/1684 


2170/1589 
3040/1710 


4587/3232 
5757/3410 


4435/3226 
6064/3404 


4818/3398 
4818/3398 



(«) w e give the separate fits for the two observations only for Model A in order to show that they are consistent within the errors and thus can be 
fitted together. C-Stat / dof (*) refers to the statistics obtained by including the wavelength ranges 17.2-17.7 A and 22.7-23.2 A, which in fits are 
ignored except in the case of the complete model. ( ' All the columns for the dabs and amoJ components are reported in units of 10 21 irT 2 . (r) Each 
amol component is displayed together with its molecular index as reported in Table [7] 
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Fig. 6. Fe i L-edge: data and Model B. 
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Fig. 7. Ne i K-edge: data and Model B. 



are described briefly below. The transmission of dust is cal- 
culated by the dabs model in SPEX: it accounts for the self- 
shielding of X-ray photons by dust grains, but uses the edge and 
line structure for atomic gas. It was first use d in the analysis 
of the Crab spectrum by iKaastra et all (l2009h. It follows com- 
pletely the dust treatment as described by Wil ms et all (l2000h 
and is useful to estimate the dust-to-gas column ratio and the de- 
pletion factor for several elements. We assume the default values 
for the grain parameters because they are physically acceptable: 
the grains are assumed to be spherical and fluffy, with density 
p = 1000 kg irT 3 , grain radius a between a m ; n < a < a max , 
where a m ; n = 0.025 fim and a m:l x = 0.25 jum, with a size distri- 
bution dn/da ~ a~ p and p — 3.5 (IKaastra et alj|200 9i). Including 
the dabs component in our model, we obtain a significant im- 



provement to the fit by requiring at least ~ 90% of iron to be 
confined in dust grains (see Table [3). Inside the iron edge the 
shielding effect of dust is stronger than its fine structure features 
and thus the dabs component is suitable to fit the data. Indeed, 
the ~ 4 cr deviation at 17.4 A (Fig. |6]) and the ~ l<x deviation 
near 17.1 A, due to the assumption of a pure-gas ISM, just dis- 
appear (see Fig. [9]). Moreover, from the dabs model we derive 
~ 40% of the oxygen to be bound in dust grains. Unfortunately, 
this model does not yield a good fit of the oxygen edge, where 
strong features due to molecules are present that are not taken 
into account by the dabs model. 

A much larger improvement is obtained when we introduce 
molecules and minerals containing oxygen atoms. For the first 
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Fig. 10. Oxygen edge: data and Model C. 
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Fig. 9. Iron edge: data and Model C. 
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time we use the amol model in SPEX to take into account 
bound forms of oxygen. The model currently takes into account 
the modified edge and line structure around the O K-edge us- 
ing measured cross sections of various compounds, taken from 
the literature. More details about this model can be found in 
Appendix [A] The amol model is very useful to constrain the 
local molecular features, but it does not account for the dust 
shielding effects. We have tried 23 different types of compounds 
such as CO, N2O, H2O, ice, FeO and several minerals. The best 
fit is obtained by using a mixture of silicates, such as andra- 
dite, and water ice plus other molecules (see Fig. [10] and Table 
|4j, This additional component removes the previous 3 <x devia- 
tion of the pure-gas model between 22.7-23.0 A inside the oxy- 
gen edge (Pig. O. Finally we complete the dust model choosing 
the amol component for oxygen and dabs component for all the 
other depleted elements, such as iron and magnesium. The final 
gas+dust model describes the data much better than all previous 
models (see Fig. [9land ITOl and it allows to estimate the dust- 
to-gas ratio. We label this model C and show the parameters in 
Table [3] 



3.3. ISM model complexity 

The spectral modeling indicates that in the line of sight towards 
our X-ray source the ISM is much more complex than simple 
neutral gas. The gas is rather structured in different phases and 
dust also consists of various compounds. 

The gas consists of 3 components (see Table [3): cold gas 
with a temperature kT ~ 5 - 10 X 10~ 4 keV (5.8-10 xlO 3 K), 
warm ionized gas with IcT ~ 1 - 5 x 10 3 keV (1-6 xlO 4 K) 
and hot ionized gas with kT ~ 0.2 keV (> 2 x 10 6 K). The col- 
umn densities A^e of these three components span over 2 orders 
of magnitude: the cold gas accounts for ~ 90 - 95 % of the to- 
tal column N l °\ N™ rm ~ 5 - 10 % of A^f, while the hot gas 
contributes ~ 1 %. 

The warm component produces the low-ionization absorp- 
tion lines of O n and O in, at 23.35 A and 23 . 1 A respectively (see 
Fig. |5j. It also provides a better modeling of the neon edge (see 
Fig. |7]i. Using model B we estimate N 0n = 2.0 + 0.5 x 10 21 irT 2 
andNom = 1.4 + 0.5xl0 21 nr 2 , respectively ~ 7% and ~ 5% of 
the total oxygen column. However, the derived column density 
of the warm ionized gas is affected by the presence of dust and 
molecules on the line of sight. Indeed, the absorption features 
that we see near 23.35 A and 23.1 A are contaminated by dust 
and molecules effects. Contributions from dust and molecules 
(Model C in Table [3} are confirmed by the improvements to the 
fit (see also Fig.l9land[T0l> and the column density of the warm 
gas is finally reduced to about 5% of the full gas column (see 
also Table |4]i . 

The column density of hot gas is about two orders of mag- 
nitude lower than the cold gas column and its temperature is 
around two million K. As expected, the hotter gas has a higher 
velocity dispersion (see Table|3]l. The hot gas model gives a good 
fit of the O vii absorption line at 21 .6 A, together with the small 
feature at 13.4 A produced by Ne ix. 

According to the analysis of the oxygen edge, the solid phase 
of the ISM towards GS 1826-238 consists of a mixture of miner- 
als (such as andradite silicates) and traces of CO and water ice. 
We cannot yet distinguish between amorphous and crystalline 
phases. As reported in Table |4] the bulk of the oxygen, ~ 90%, 
appears to be in the gas phase, while the remaining ~ 10% is 
made mostly of solids, such as silicates and water ice. Obviously, 
there could be substances able to reproduce such features in the 
spectrum other than our few dozen test molecules. For the iron, 
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Table 4. 

density. 



Table of the contributions to the oxygen column- Table 6. Ionic column densities of oxygen in 10 22 m 2 . 



Phase 


Constituent 


No (1CF m" 2 ) 


% of JV£ <"' 


% of N {b) 


Gas 


Oi 

On, m, O iv 

O vii, O vm 


2.7 ±0.1 
0.10 ±0.05 
0.05 ±0.01 


94 ±4 
4±2 
2.0 ±0.5 


90 ±6 


Dust 


Silicates 
Aluminates 


0.25 ± 0.05 
<0.04 


85 - 100 
0-15 


10 ±2 


Molecules 


H 2 ice 
CO 


<0.07 
<0.04 


-65 

-35 


0-2 



(a) % of Nq represents the contribution of each constituent to the respec- 
tive phase. (6) % of A'o give the contribution of the different phases to 
the total oxygen column density. See also Sect. l3,2,3l and Fig.[TT|for the 
transmission of the main compounds. 
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Fig. 11. Transmission near the Oi K-edge. The model used is 
Mod C (see Sect. 1X231 Tableland©. Both cold-gas and entire- 
ISM transmission is multiplied by a factor of 4 for displaying 
purpose. 



instead, we obtain a different composition: at least ~ 90% of Fe 
appears to be bound in dust grains. In our dust model we as- 
sume a depletion factor of 0.8 for magnesium, as suggested by 
Wil ms et alJ (|2000). The derived gas (2.0 + 0.4 x 10 21 trT 2 ) and 
dust (2.3 + 0.1 x 10 21 irT 2 ) column densities for the Mg are iden- 
tical within the errors (see also Table[3]). 



3.4. ISM abundances 

We have estimated the abundances of Mg, Ne, Fe, O and N 
(Table |5). The column density of each element refers to the sum 
of the contributions from all the gas and dust phases. The abun- 
dances do not differ significantly between the two observations 
(see TableO. This is expected if the absorption is mainly due to 
the interstellar medium, because the ISM is stable on short time 
scales. The zero shift of the Oi line, the position of the O, Fe, 
Ne and Mg edges (see Fig. 151 to [Tot and the low velocity dis- 
persion suggest that the absorber matter is a mixture of gas and 
dust without outflows or inflows, not broadened due to Keplerian 
motion around the X-ray source. This is also consistent with an 
ISM origin. For further discussion on the abundances and com- 
parisons with previous work see Sect. 14.31 



X 


Mod A,„„, (fl| 


Mod A (fl) 


Mod B m 


Mod C (r) 


Oi 

O ii, O in, O iv 
O vii, O vm 
Dust 


3.0 
= 
= 
= 


3.1 

= 
= 
= 


2.9 
0.35 
0.05 
= 


2.7 
0.1 
0.05 
0.35 


Total A'o 


3.0 


3.1 


3.3 


3.2 



<fl) Mono-phase gas (Table [3] for the simultaneous EPIC-RGS fit). (i) 
Three-phases gas (Table |3j- (c) Gas + dust model (see Table |4j- The 
agreement between the oxygen column densities estimated with differ- 
ent models validates our method. 



A separate analysis is required for the nitrogen edge at 30.8 
A. Indeed, at wavelengths higher than 29 A the source flux de- 
creases significantly and is comparable to the background. Thus, 
the spectrum around the N i K-edge is noisy and the N i column 
density has a large uncertainty. We have decided to freeze the 
continuum parameters and try either a broadband fit to the range 
7-33 A or a local fit between 27-33 A. The abundance estimates 
are consistent within the error and their average is reported in 
Table|5j 

4. Discussion 

4.1. The continuum 

Our analysis shows that the persistent state of the low-mass X- 
ray binary GS 1826-238 is well represented by a double comp- 
tonization (C1+C2) plus a black body (BB) emission compo- 
nent, all three absorbed by the interstellar medium composed of 
a 3-phases gas, dust and molecules (see Table [2] and [3}. In both 
observations the fits are in agreement: all the ISM parameters ap- 
pear to be fully consistent, thus we can discuss about the results 
from the simultaneous fit of the high-resolution RGS data. 

4.2. ISM structure 



In Sect. l3.3l we show that in our line of sight the ISM has a clear 
multiphase structure. There are media with different ionization 
states, dust grains and molecules. As confirmed by Fig. [TT] the 
bulk of the matter responsible for X-ray absorption is found in 
the form of cold gas with a temperature - 7 000 K and low ve- 
locity dispersion (cry < 13 km s -1 ). At this temperature the gas 
is almost neutral: only iron and magnesium are partially ionized. 
Part of the cold matter is found in solids compounds, such as 
dust grains and molecules. Most of iron is bound in dust grains. 
About 10 % of oxygen is found in compounds: the silicates con- 
tribute to ~ 80% of this phase, while the remaining fraction 
consists of a mixture of other oxides (such as iron aluminates) 
together with ices and CO molecules in similar quantities (see 
Table @J. The best fit is obtained using as compound the andra- 
dite silicate Ca3 Fe2 (SiOi)3, but we need higher signal-to-noise 
data to distinguish among the different silicates, as also olivine 
and pyroxene are good candidates. Moreover, at the present stage 
our model does not take into account simultaneously the shield- 
ing and fine structure effects of oxygen compounds, thus the 
fraction of oxygen bound in solids could be higher, e.g. up to 
40 % (see Sect. 13.2.3b . However we are working to the devel- 
opment of models that take into account all the possible effects 
and we postpone a deeper analysis of the oxygen dust phase to a 
forthcoming paper. 
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Table 5. Average ISM abundances in units of the proto-Solar values (Lodders 2003) and calculated by summing contribution from 
all ISM phases. 



X 


O 


Ne 


Mg 


Fe 


N W 


GS 1826-238 {a) 


1.23+0.05 


1.75+0.11 


2.45 ± 0.35 


1.37 ±0.17 


2.4 ± 0.7 


Crab (A) 


1.030 ±0.016 


1.72 ±0.11 


0.85 ± 0.21 


0.78 ± 0.05 


1.01 ±0.09 


Cyg X-2 (c) 


0.6-0.8 


0.8- 1.1 


0.6-1.1 


- 


- 


4U 1820-303 w 


0.7-1.1 


1.1-2.0 


- 


0.3 - 0.8 


- 



<"' Model C, gas + dust. m lKaastra etafl d2009h Model B, gas + dust. (c) lYao et alJ J2009l> . w Uuettet all 12006). 
{e) Estimated through a local fit in the 27-33 A range. 



About 5% of the gas is ionized (see Mod C in Table[3]l. Most 
is a warm plasma with T ~ 10-50 000 K. It has a low-ionization 
degree and it accounts for the interstellar On and Om. Only 
1-2% of the ISM gas appears to be highly ionized. Such a hot 
plasma reaches temperatures ~ 2 x 10 6 K and contributes to all 
O vn and O vm present in the ISM. As we expect, the higher the 
temperature of the gas phase, the higher its velocity dispersion. 
Unfortunately the velocity dispersion is not well constrained, es- 
pecially for the hot ionized gas. This is not surprising as the lines 
are unresolved. 

4.3. Comparison with previous results 

4.3.1. ISM constituents 

The average total column density of the multi-phase gas we esti- 
mate is about(4.14 + 0.07)xl0 25 irT 2 . This is not consistent with 
the va lue of (3.19 + 0.01) x 10 25 irT 2 found bv lThompsonet all 
(l2008h . They combined data from Chandra, XMM-Newton and 
RXTE in the 0.5 - 100 keV band, but we know that the bulk 
of the absorption is at lower energy. Instead we have used both 
EPIC and RGS data. The latter detector has higher spectral 
resolution between 0.3 - 2 k eV and allows to better estimate 
the absorption column. Also, Thompson et all (2008) used the 
Morrison & McCa mmonl (I 1 9831) model to fit absorbing medium, 
which takes into account only contribution by cold neutral gas. 
In Sect. I3.1l we have shown that such a model is not optimal, be- 
cause it does not include lines, then abundances and temperature 
cannot be free parameters. Instead our estimate of Nh is obtained 
by summing contribution from all the phases of the gas and by 
accounting for all the absorption features found in the spectrum. 

The multiphase structure of the ISM that we have con- 
strained is consistent with rece nt results (iFerrierd 1200 it 
lYao & WangH2006t lYao et alj|2003) . First of all, there is a good 
agreement in the fractions between the cold, warm and hot 
phases of the gas. In particular the estimated amount of O vn, 
~ 1.6 x 10 20 m~ 2 , is fully consistent with the value found by 
lYao & Wangl (120061) by fitting both O vn ls-2p and ls-3p lines 
in the spectrum of 4U 1820-303, which is another LMXB near 
the center of the Galaxy. The hot gas accounts for 1 - 2 % of 
the total column density and it represents the average fraction 
of hot plasma in the Galaxy. Indeed it agrees with previous es- 
timates and can be fully expla ined by the heating of supernovae 
(McCammon & Sanders 1990). The velocity dispersion cry es- 
timates agree with recent work. Yao & Wang (2006) combined 
oxygen and neon ionization lines in the Chandra spect rum of the 
LMXB 4U 1820-303 obtaining cr v < 350 km s" 1 . Juett et al. 
(2004) found cry < 200 km s _1 from oxygen lines fits to the 
Chandra spectra of several LMXBs. 

Moreover, we found clear indications of dust depletion in 
some heavy elements, such as oxygen, iron and magnesium. 
According to the different dust models that we have used (dabs 
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Fig. 12. Map of the X-ray sources compared in this paper. GC 
is the Galactic Center and the Sun is assumed to be 8.5 kpc far 
away from it. 



and dabs+amol), we find ~ 50% of Mg, more than 90% of Fe 
and 10 - 40% of oxygen in the form of dust grains and molecu- 
lar compo unds. Such results are mostly consisten t with previous 
estimates dWilms et al.l2 000: Kaastra et al. 2009). The iron dust- 
to-gas ratio found towards GS 1826-238 is among the highest 
measured in the Galaxy. Williams & Taylor] d 19961) showed that 
the higher is the density of a molecular cloud, the higher is the 
probability of forming dust and molecules from gas particles. 
Thus the higher dust-to-gas ratio we estimated towards the cen- 
ter of the Galaxy suggests high density regions. 

The presence of silicates and ice in the ISM, c onstrained by 
the O i K-edge analysis, is supported by other work. lPaerels et al] 
(2001) found similar features within 22.7 - 23.0 A in the spec- 
trum of the LMXB X06 14+091 and they argued that it should 
be due to iron oxide s or oxygen generally bound in dust. 
ICostantini et al.l d2005l) found indications of silicates such as 
olivine and pyroxene by modeling the feature of the scattering 
halo of Cyg X-2. Recently. Ide Vries & Cost antini (2009) found 
evidence for EXAFS in the short-wavelength side of the oxygen 
edge in the spectrum of Sco X- 1 and their results suggest the 
presence of amorphous water ice. 

We have checked CO surveys (Da me et aD l2001) in order to 
test the presence of molecular clouds. There is no clear evidence 
for these clouds in our line of sight (x,y)cAL ~ (9.3°, -6.1°), 
while the re is an im portant amount of diffuse dust (see also 
Schlegel et al. 1998). This is consistent with the fact that we can 
only put upper limits to ice and CO columns. Thus the absorption 
should arise from a solid phase consisting mostly of minerals. 
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4.3.2. ISM abundances 

The abundances are displayed in Table [5] they are reported rela - 
tive to the recommended proto-Solar values of Lodders (2003). 
We also show the abundances obtained by Juett et al. (2006) to- 
wards the LMXB 4U 1820-303, which like GS 1826-238 is 
close to the Galactic center. Even more int eresting is the com- 
parison with the abundances estimated bv iKaastra et al. (2009) 
in the direction of the Crab nebula and the ones estimated by 
Yao et al. (2009) towards Cyg X-2, which are two different lines 
of sight in our Galaxy (see the map in Fig.[T2l. 

The derived abundances slightly depend on the model used: 
the pure-gas model and the complete gas+dust model yield dif- 
ferent results on Ne, Fe and Mg (Table[3]l. First, their abundances 
are more uncertain than the oxygen abundance. Secondly, Fe and 
Mg are the most depleted elements: as dust grains give rise to ab- 
sorption features different from gas particles, modeling edges of 
highly depleted elements provides different results by account- 
ing or not for dust effects. Such deviations are probably strength- 
ened since we observe throu gh a high-density region : here high 
metal depletion is expected (Williams & Tavlorl fl996). 

All the heavy elements are over-abundant with respect to the 
proto-Solar values. Neon is over-abundant by a factor ~ 1.7, as 
found by Kaastra et al. (2009) in the XMM-Newton observations 
of the Crab nebula. The Solar abundance of neon is probably 
under-estimated (see lLodders & Palmd 120091) . so our estimate 
may be not really different from the real Solar value. As we will 
show in the next paragraph, the metallicity gradient could also 
be responsible for part of the Ne over-abundance. 

The reason for the over-abundances of O, Fe and Mg is rather 
different from the neon excess. First of all, we are able to mea- 
sure both gas and dust contributions for O, Fe and Mg, without 
the risk of missing important fractions. Secondly, deviations in 
the abundances with respect to the average Galactic values are 
also due to their metallicity gradient. If A(X) is the abundance 
of a certa in element X in the vicinity of GS 1826-238, we can 
write (see Esteban et al. 2005), 



A(X) 
A(X) Q 



= 10 



i«x(Ogs-D ) 



(1) 



where D Q and Dqs are the Galactocentric radii of the Sun and 
GS 1826-238, respectively ~ 8.5 kpc and ~ 2 kpc; ax is the 
abundance gradient of the element along the line of sight and 
A{X) G is its abundance near the Sun. In this way, we can compare 
our estimates of abundance changes with the values predicted 
by the gradient estimates. Unfortunately, the gradient estimates 
in the literature are quite uncertain and are availabl e only in a 
limited range of radii, i.e. between 4-16 kpc (see Pedicelli et al. 
2009). Thus we can trust only in the abundance changes in 
4-5 kpc along our line of sight. Moreover, the column den- 
sity estimated for each element refers to its integral along the 
line of sight, where we also expect a density increase towards 
the Galactic Center. Because the density increases towards the 
Galactic center, the predicted abundance at the Galactocentric 
distance of GS 1826-238 should be close to the weighted aver- 
age abundance along the sightline: 

- Oxygen is over- abun dant by about 20 - 30 % (see Table[5]). 
According to lEsteban et al.l (|2005), the oxygen gradient is 
ao - (-0.04 + 0.01) kpc -1 , which should provide an incre- 
ment (equation[T} of at least ~ 32% in the oxygen abundance. 
This is consistent with our estimate. 

- The iron abundance ~ 1.20-1.54 (see Tabled is near agree- 
ment with the > 50% increment derived by the accepted iron 



gradient in the Galactic disk ap e = (-0.06 ± 0.02) kpc 
(iFriel et al.ll2002t iPedicelli et al.ll2009h . 

- Neon is over-abundant by more than 70%. It is difficult 
to compare it with the Galactic gradient because this is 
quite uncertain, on average q-Ne = (-0-06 + 0.04) kpc ~' 
(Simp son et alj|1995t iMaciel & Ouirezalll999l) . From equa- 
tion ([TJ we predict a lower limit A(Ne) ~ 1.3. The sum 
of such a value to the revisited proto-Solar abundance 
(AA(Ne) ~ 30 %, see lLodders & Palmd [2009b provides 
A(Ne) ~ 1.6, which is fully consistent with our estimate. 
This result suggests that the neon over-abundance that we 
constrain is due to both the Galactic gradient and the previ- 
ous under-estimate found in the literature. 

- Nitrogen show s a steeper Galac t ic gra dient of about ~ 
-0.08 kpc 1 (Gummers bachet all 119981) . which provides 
A A(N) ~ 100 % (equation [TJ and agrees with our estimate 
(see Table[5]l. Of course, both of them are quite uncertain and 
we are not able to provide more information. 

- Also magnesium should be at least twice the proto-Solar 
value as we found (see Table [5}. However, it is difficult to 
compare our result with the Mg Galactic gradient found in 
the literature, because the results differ a lot in the literature 
(Roll eston et al.l 120001) . According to Gummersbac h et al.l 
( 1998) a Mg ~ -0.08 kpc 1 , which together with equation CO 
implies abundance increment of at least 100%. This is fully 
consistent with what we found. 

The iron and oxygen abundances towards 4U 1820-303 
(Table [5) appear to disagree with our results des pite its similar 
location near the Galactic center. Hue tt et al.l (2006) attribute this 
low iron abundance to depletion into dust grains in the interstel- 
lar medium and they also report that oxygen could be middle- 
depleted. In Sect. 13.2.31 we have shown that iron is among the 
most depleted elements in the ISM, oxygen is also partially de- 
pleted and a pure-gas model cannot reproduce all the ISM spec- 
tral features. This indicates that the differences between the Fe 
and O abundances are due to the capability of our gas+dust 
model to measure the contribution from solid phases, which are 
absent in the pure-gas model of I Juett eta l. (2004, 2006). 

The Crab nebula is relatively close to the Solar system, i.e. 
~ 2 kpc, but opposite to the Galactic center (see Fig. ITZI ). Thus 
we expect abundances clo se to the proto-Solar v alues of lLoddersI 
(2003), which is just what lKaastra et alj (120091) found. 

LMXB Cyg X-2 is also far away from the center of the 
Galaxy and about 10 kpc away from the Sun. The abundances 
estimated for this source are lower than those measured towards 
the Crab and GS 1826-238 (see Table [5). This agrees with the 
assumed abundance gradients we have discussed. 

In summary, the differences with respect the proto-Solar 
abundances that we estimate are consistent with the literature. 
The increase of metallicity towards the center of the Galaxy 
should be due to evolutionary effects like supernovae explosions, 
which enrich the ISM with heavy elements, especially in the 
higher density region of the bulge and the disk of our Galaxy. 
A deeper analysis is required: we need to study more sources, 
even in the same region, and further improve our models in or- 
der to account for every contribution to the column densities. 



5. Conclusion 

We have presented a complete treatment of the interstellar 
medium towards the low-mass X-ray binary GS 1826-238, 
which is a bright X-ray source near the Galactic center. We have 
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shown that in the line of sight the ISM is composed of a complex 
mixture of a multi-phase gas, dust and molecules. 

The gas is almost neutral and the ionization degree is about 
5%. Significant fractions of the column density of some heavy 
elements are in the form of molecules or dust grains: at least 
10% of oxygen, 50% of magnesium and 90% of iron. Such a 
solid phase should consist of a mixture of silicates (> 60%), CO 
(< 10 %), ice (< 20 %) and other iron oxides (< 10 %). 

We have found over-abundances for all the elements for 
which we have been able to measure the individual column den- 
sity. The Ne over-abundance that we estimate is consistent with 
other rec ent work in d ifferent lines of sight, such as towards 
the Crab dKaastra et alj |2009). suggesting that the Solar value is 
underestimated. However, our agreement with the predicted Ne 
gradient in the Galaxy is also consistent with an abundance in- 
crease due to stellar evolution: towards the Galactic center there 
are high-density regions with evolved star that could have en- 
riched the ISM with heavy elements such as neon. 

Differently from the previous X-ray spectroscopy work, we 
have found over-abundances for oxygen (1.2), iron (1.4) and 
magnesium (2.4). These elements are also in the form of dust 
and molecules. Thus our estimates are partially due to the fact 
that we are able to measure also the contributions from the solid 
phase. The abundance of metals is in agreement with the metal- 
licity gradients and shows the chemical inhomogeneity of the 
interstellar medium. 

The diagnostic of the ISM constituents fits the predicted 
models for both its thermal and chemical structures, showing a 
good agreement with the current state of art (Ferriere 2001). The 
dust column is consistent with the multi-wavelength, X-ray ver- 
sus IR, observations. All this supports our research method and 
justifies new observations and analysis of other sources in dif- 
ferent lines of sight. Such analysis, indeed, can provide a better 
mapping of the ISM and a deeper study of its chemical composi- 
tion, together with its role in the evolution of the entire Galaxy. 
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A. Amok oxygen edge molecules absorption model 

The amol model calculates the transmission of various 
molecules. Presently only the oxygen edge is taken account of. 
Updates of this model, once made, will be reported in the man- 
ual of SPEX. The following compounds are presently taken into 
account (see Table|7]i. 

The chemical composition of these minerals was mainly 
taken from the Mineralogy Database of David BarthelmyQ We 
take the cross-sections from the references as listed in TableQin 
the energy interval where these are given, and use the cross sec- 
tion for free ato ms Verner & Yakovlev ( 1995) outside this range. 
Ivan Aken et alj (1 19981) do not list the precise composition of iron 
oxide. We assume here that x - 0.5. 

Some remarks about the data from lBarrus et al.l (1 19791) : not 
all lines are given in their tables, because they suffered from in- 
strumental effects (finite thickness absorber combined with fi- 
nite spectral resolution). However, Barms et al. have estimated 
the peak intensities of the lines based on measurements with dif- 
ferent column densities, and they also list the FWHM of these 
transitions. We have included these lines in the table of cross 
sections and joined smoothly with the tabulated values. 

For N2O, the fine structure lines are not well resolved by 
Barms et al. Instead we take here the relative peaks from 
Wig ht & Brionl (1 19741) . that have a relative ratio of 1.00 : 0.23 
: 0.38 : 0.15 for peaks 1, 2, 3, and 4, respectively. We adopted 
equal FWHMs of 1.2 eV for these lines, as measured typically 
for line 1 from the plot of Wight & Brion (1974). We scale the 
intensities to the peak listed by B arms et aD (Il979l) . Further, we 
subtract the C and N parts of the cross sectio n as well as the oxy- 
gen 2 s/2p part, using the cross sections of 1 Verner & Yakovlevl 
(1995). At low energy, a very small residual remains, that we 
corrected for by subtracting a constant fitted to the 510-520 eV 
range of the residuals. The remaining cross section at 600 eV is 
about 10 % above the Verner cross section; it rapidly decreases; 
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Table 7. Molecules present in the amol model. 



chemical formula reference 



1 


molecular oxygen 


o 2 


a 


2 


carbon monoxide 


CO 


a 


3 


carbon dioxide 


co 2 


a 


4 


laughing gas 


N,0 


a,b 


5 


water 


H,0 


c 


6 


crystalline ice 


H,0 


d 


7 


amorphous ice 


H 2 


d 


8 


cupric oxide 


CuO 


c 


9 


nickel monoxide 


NiO 


c 


10 


iron oxide 


Fe,_. r O 


e 


11 


magnetite 


Fe 3 4 


e 


12 


hematite 


Fe 2 3 


c 


13 


eskolaite 


Cr 2 3 


c 


14 


andradite 


Ca 3 Fe 2 Si 3 12 


c 


15 


acmite 


NaFeSi 2 6 


c 


16 


franklinite 


Zno.6Mno.gFe1.6O4 


c 


17 


chromite 


FeCr,0 4 


e 


18 


ilmenite 


FeTi6 3 


e 


19 


perovskite 


CaTi0 3 


c 


20 


olivine 


Mg 1 .6Feo. 4 Si0 4 


c 


21 


almandine 


Fe 3 Al 2 (Si0 4 ) 3 


c 


22 


hedenbergite 


CaFeSi 2 06 


c 


23 


hercynite 


FeAl 2 4 


e 



r eferences: 

Barrus et al.Ul979l). .5-0.75 eV resolution 
Wight & Brionl <1974h . 0.5 eV resolution 
Hirava et al.1 (l200lh . 0.055 eV resolution 
Parent et al.fd2002h. , 1 eV resolution 
van Akenetal.Ul998). 0.8 eV resolution 



we approximate the high-E behavior by extrapolating linearly 
the average slope of the ratio between 580 and 600 eV to the 
point where it becomes 1 . 



